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It is shown that spatially periodic one-dimensional surface waves in shallow water behave almost lin- 
early, provided large part of the energy is contained in sufficiently high frequencies. The amplitude is not 
required to be small (apart from the shallow water approximation assumption) and the near-linear behavior 
occurs on a much longer time scale than might be anticipated based on the amplitude size. Heuristically 
speaking, this effect is due to the nonlinearity getting averaged by the dispersive action. 

This result is obtained by an averaging procedure, which is briefly outlined, and is also confirmed by 
numerical simulations. 



The dynamics of surface water waves has been an im- 
portant object of study in science for over a century. Soli- 
ton solutions and integrability in P.D.E.'s are two examples 
of remarkable discoveries that were made by investigating 
water wave dynamics in shallow waters. In more recent 
times, the so-called rogue waves have been under an in- 
tense investigation, see for example HI d, 01 and the ref- 
erences therein. These unusually large waves have been 
observed in various parts of the ocean in both deep, see 
e.g. Jl, and shallow water, see e.g. M, 

leading scientists 

to suggest various mechanisms for rogue wave formation. 

In the case of shallow water, one normally does not 
work with the full water wave equation but uses approx- 
imate models to study the evolution, in particular rogue 
waves. These models are nonlinear dispersive equations 
such as KdV, Boussinesq approximations, etc. In partic- 
ular, KdV describes unidirectional small amplitude long 
waves on fluid surface. See, e.g. [8] for applications of 
KdV to rouge waves in shallow water. Since rouge waves 
correspond to concentration of energy on small domains, 
it is expected that high frequencies play important role in 
rouge waves formation. 

In this Letter, we demonstrate that for sufficiently high 
frequency initial data, one-dimensional spatially periodic 
surface waves in shallow water exhibit near-linear behav- 
ior. This phenomenon is established by observing that high 
frequency periodic solutions in KdV equation behave lin- 
early. 

Thus, linear theories of rogue wave formations can be 
extended to nonlinear high frequency regime. While, we 
demonstrate this phenomenon for KdV, we expect that sim- 
ilar behavior will occur for other models describing shal- 
low water waves. We note that our near-linear dynamics 
has nothing to do with weakly nonlinear theories as we do 
not require small amplitude or small energy. The evolution 
is near-linear on larger time-scale than one might expect 
based on the size of nonlinearity. Moreover, we deal with 
dynamics on bounded domains (periodic boundary condi- 
tions) so the solutions cannot scatter to infinity. Near-linear 
behavior for high frequency solutions on unbounded do- 
mains is easier to establish as, essentially, the energy dis- 
perses to infinity and weakly nonlinear theories become ap- 
plicable. 



By contrast, on bounded domains {e.g. periodic, Dirich- 
let, or Neumann boundary conditions), near-linear be- 
havior may occur only if the net effect of nonlinearity 
somehow averages out, as the nonlinearity stays of finite 
strength. Heuristically, one may speculate that strong dis- 
persion will cause high harmonics to oscillate rapidly and 
average the nonlinearity out, assuming that there is no fo- 
cusing collapse in the problem. While such explanation 
seems plausible, it is important to formulate precisely un- 
der what conditions and on what time scales such near- 
linear dynamics takes place and this is what we undertake 
in this work. 

Clearly, one has to be careful when considering short 
wave solutions for the equations obtained in the long wave 
approximations such as KdV. However, we show that there 
is a set of parameters when our high frequency solutions 
correspond to a realistic physical scenario in shallow water 
waves, see Section. 



NEAR-LINEAR DYNAMICS AND ROGUE WAVES 

Rogue waves are large-amplitude waves appearing on 
the sea surface seemingly "from nowhere". Such abnor- 
mal waves have been also observed in shallow water and 
KdV has been used to explain this phenomenon [8j]. In 
the oceanographic literature, the following amplitude cri- 
terion for the rogue wave formation is generally used: the 
height of the rogue wave should exceed the significant 
wave height by a factor of 2-2.2 01. (Significant wave 
height is the average wave height of the one-third largest 
waves.) 

There is a vast literature on rogue (also called freak and 
giant) waves, as discussed above, and many explanations 
have been proposed. Major scenarios and explanations in- 
volve 

• probabilistic approach: rogue waves are considered 
as rare events in the framework of Rayleigh statistics 

• linear mechanism: dispersion enhancement (spatio- 
temporal focusing) 

• nonlinear mechanisms: in approximate models {e.g. 
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KdV), for some special initial data, large amplitude 
waves can be created. 

Linear mechanism is very attractive as there are sim- 
ple solutions leading to large amplitudes, while nonlinear 
mechanism requires rather special initial data, e.g. leading 
to the soliton formation. On the other hand, linear equa- 
tions arise in the small amplitude limit which is too restric- 
tive. 

Using near-linear dynamics in KdV one can experiment 
with a new mechanism of large wave formation that com- 
bines linear and nonlinear deterministic mechanisms. Our 
results indicate that for a special but relatively large set of 
initial data (characterized by the energy contained mostly 
in high frequency Fourier modes), the solutions of KdV 
equation behave near-linearly. It is then possible to con- 
struct large amplitude solutions using linear mechanisms 
of large wave formation. 

Here, we illustrate our approach with periodic boundary 
conditions. This is not the most realistic choice but appro- 
priate for a model problem to illustrate the concept. Indeed, 
while the sea surface is not periodic, one observes similar 
patterns over large areas. It should be possible to extend 
near-linear dynamics to quasiperiodic and random initial 
data and we intend to do it in the future investigation. 

It would be also desirable to observe near-linear dynam- 
ics for the full water wave problem, however, it is a con- 
siderably harder problem which we also hope to address in 
our future work. 

The fact that nonlinear dispersive systems exhibit near- 
linear behavior seems to be rather general and should be 
important in other fields where nonlinear dispersive equa- 
tions appear as various approximations, such as solid state 
physics, photonics, etc.. 



MECHANISM OF NEAR-LINEAR DYNAMICS IN 
NONLINEAR DISPERSIVE SYSTEMS 

Nonlinear dispersive equations that are derived as short 
or long wave, weakly nonlinear approximations of basic 
physical models, often take the form 

u t = Lu + f(u) 

for a field variable u which may represent slowly varying 
amplitude, rescaled velocity, etc. The operator L is a linear 
differential operator and / is nonlinearity which may also 
contain derivatives but of lower order than those in L. 

On infinite domains, dispersion tends to spread localized 
waves whereas nonlinearity may cause shocks or focusing. 
Here, as discussed in the introduction, we consider the phe- 
nomenon of dispersive averaging when the dynamics oc- 
curs in a finite region and dispersion cannot cause energy to 
spread to infinity. It turns out that dispersion creates rapid 
oscillations for high harmonics which effectively average 
out nonlinearity. If most of the energy is contained in high 



frequencies, then Lu will dominate over f(u) assuming 
that we consider subcritical nonlinearity. However, it is not 
clear what is the effect of the nonlinearity over large times 
since u is not small. We show that the nonlinear effects are 
small leading to essentially linear evolution for finite but 
large times. 
Let 



u 



e Lt v, 



so that 



v t = e- Lt f(e Lt v). 



If e is periodic in time then one can attempt to apply the 
usual averaging procedure directly. However, it is not clear 
at all that the averaging procedure applies in this setting. 
Indeed, the above system is a partial differential equation 
and also there is no clear separation of time scales. Nev- 
ertheless, our main result described in the next section and 
subsequent numerical simulations indicate that it could be 
possible to get around these difficulties and show that the 
above system is close to the averaged one 

v t = (e- LT f(e LT v)) T . 

Now, we describe this mechanism for KdV but, as we have 
mentioned, it seems to be applicable with appropriate mod- 
ifications to many other systems. 



NEAR-LINEAR DYNAMICS IN KDV 



Consider KdV 



u t = 6uu x + u xxx 



with periodic boundary conditions n(x + 27r) = u{x). The 
following statement characterizes near-linear dynamics in 
this equation: 



Assume that initial data u(x) has zero mean (u) = 
and has energy equal to 1: 



/u 2 (x)dx = 
-7T 



and assume that Fourier transform u satisfies 



E 

n=£0 



u(n) 



<Ce 2 . 



Then, 



\u(x,t) - e Lt u(x,0)\ 2 dxj < C 2 (e 2 + te), 

where e Lt is the free propagator for KdV. 

The proof of this statement will appear elsewhere and 
here we only present heuristic arguments. 
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In Fourier domain, 

u (x) =Y J v{n)e mx 

n 

the equations take the form 

v(n) = in 3 v(n) + 6 n 1 w(n 1 )t>(n 2 ). 

ni+n 2 =n 

After a time dependent change of variable 

v(n) = e m3t w(n) 

we obtain 

w(n) =6 e min2nt n lW {ni)w{n 2 ). 

n 1 +ri2=n 

The goal is to show that w stays almost constant for large 
times under our high frequency assumptions on the initial 
data. It is easy to see that for general solutions, one cannot 
expect this kind of behavior and on the other hand, we can- 
not estimate rate of change of w by neglecting averaging 
effects of the exponent. Indeed, without the exponential 
factor e miIl2nf , the above system corresponds to Burgers 
equation u t = uu x which is known to have shocks and is 
known to exhibit strongly nonlinear dynamics. 

Using high frequency assumptions and averaging effect 
of the exponent, we find that near-linear dynamics can be 
established by applying a version of the normal form pro- 
cedure, see e.g. [5]. The detailed analysis is rather subtle 
and will appear elsewhere. Instead, we will give an outline 
of the analysis and we will present results of our numerical 
simulations. 

There are three types of terms in the above equations: 

1. Low harmonics with ni,n 2 small. They give small 
contribution because w is small in a certain norm and there 
are few such terms. 

2. High frequency harmonics with ni,n 2 large. They 
are well averaged by the exponent. 

3. Intermediate terms. There are not so many interme- 
diate terms as the Diophantine equation n\n 2 n = N has 
few solutions ( on the order of log N). 

As a result, when one carries out two normal form trans- 
formations, which correspond to changing variable w so 
as to absorb the action of non-resonant terms, one obtains 
the system that indeed undergoes small change during the 
evolution. 



NUMERICAL SIMULATIONS AND PHYSICAL 
INTERPRETATION 

The KdV equation has been used to describe surface 
water waves in the small amplitude limit of long waves 
in shallow water. More precisely, two parameters are as- 
sumed to be small and equal 



amplitude / depth 
depth \ wavelength 




Therefore, one must be careful when considering high fre- 
quency limit in KdV as it may have little relation to the 
actual wave dynamics. However, in our case the near- 
linear dynamics phenomenon occurs when small parameter 
e characterizing high frequency limit (see the formula be- 
low), is only moderately small e = 0.4, as numerical sim- 
ulation show (one observes near-linear dynamics for finite 
time). 

As initial data we take the scaled 1st Hermit function 




so that the energy J u 2 dx does not depend on e and is very 
close to 1. For numerical simulations, we use KdV in the 
form 

3 1 

L/<£ 'XJb'XJb r£ I 'XJj 

2 6 

as it appears in the derivation of KdV from the water wave 
equations (see below). Specific numerical parameters are: 
the length of periodic domain L = 2tt. The number of 
modes M = 2 9 . Time step size At = with the time 
of the evolution T = 1. The discretization in space is given 
by h = L/M. We used the so-called Fornberg-Whitham 
scheme which is described in 111 Oh . 

Next, using standard derivation of KdV from water 
waves equations, we recall the relation between physical 
parameters and rescaled dimensionless variables, see 19fl, 
Chapter 13.11. 

Let h be the depth when the water is at rest and let 
Y = ho + rj be the free surface of the water. Let a be 
a characteristic amplitude and I be a characteristic wave 
length. Assume that 



Both a and /3 are small parameters in the problem and they 
must be of the same order. 

Next, use the following natural normalization 

x' = Ix, Y' = h Y, t' = lt/c , rf = ar], 




FIG. 1 . Initial and evolved waves in KdV in Fourier space (left) 
and physical space (right). Blue curve represents initial data 
and green curve represents nonlinearly evolved data with sub- 
sequent backward in time linear evolution, e~ Lt u(x, t). Abscissa 
shows the number of Fourier harmonic (left) and spatial coordi- 
nate (right), while ordinate is the amplitude. The time of evolu- 
tion is T = 1. 
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where primed variables are the original ones and c = 

The formal asymptotic expansion leads to KdV with 
higher order corrections 

3 1 

m + Vx + o a * + cP^xxx + 0{a 2 + /3 2 ) = 0. 
2 o 

Let X = x — t and T = at, so the equation becomes 
3 1 

Vt + TiVVx + t:Vxxx + 0(a + /3 2 /a) = 0. (1) 
2 b 

One should expect that this approximation has accuracy of 
the order 0(a) for finite time T = 0(1), which implies 
t ~ a^ 1 and t' ~ //(c a). 

Finally, since we modify our solution with another pa- 
rameter e, we verify that KdV approximation will still 
make sense for some choice of the parameters. 

First, let a = (3 = 5 -C 1. Let us modify a and I with 
a e = -^a and l e = el which is consistent with our scaling 
of initial data. Then, we have 

_ a* S 

These are small with 5 = 0.01 and e = 0.4. On the other 
hand the "mismatch" in the equation ([TJ is 

ae + ^~^ + ^~4- 

Therefore, our high frequency regime may approximate 
water waves dynamics for example with the following 
parameters: a = 1 m, ho = 100 m, and I = 1000 m. 

In summary, we have observed a new mechanism 
of near-linear behavior in strongly nonlinear dispersive 
equations. Precise formulations are given for KdV and 
implications for water waves dynamics is found. 
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FIG. 2. This figure shows the initial data and the solution after 
t = 0.2. Because the initial data is spectrally broad, the initial 
wave disperses over the whole periodic domain. 



